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Abstract. We propose to make use of the off-shell recursive relations with the color-flow decomposition in 
the calculation of QCD amplitudes on MadGraph. We introduce colored quarks and their interactions with 
nine gluons in the color-flow basis plus an Abelian gluon on MadGraph, such that it generates helicity 
amplitudes in the color-flow basis with off-shell recursive formulae for multi-gluon sub-amplitudes. We 
demonstrate calculations of up to 5-jet processes such as gg —¥ 5g, uu — > 5g and uu — > uuggg. Although 
our demonstration is limited, it paves the way to evaluate amplitudes with more quark lines and gluons 
with Madgraph. 



1 Introduction 

The Large Hadron Collider (LHC) has been steadily run- 
ning at 7 TeV since the end of March/2010. Since LHC 
is a hadron collider, understanding of QCD radiation pro- 
cesses is essential for the success of the experiments. To 
evaluate new physics signals and Standard Model (SM) 
backgrounds, we should resort to simulation programs which 
generate events with many hadrons and investigate ob- 
servable distributions. In each simulation, one calculates 
matrix elements of hard processes with quarks and gluons, 
and this task is usually done by an automatic amplitude 
' calculator. 

MadGraph[T] is one of those programs. Although it is 
highly developed and has ample user-friendly utilities such 
\ as event generation with new physics signals matched to 
Parton Shower 2 , it cannot evaluate matrix elements with 
more than five jets in the final stat43[3 . 

This is a serious drawback since exact evaluation of 
multi-jet matrix elements is often required to study sub- 
structure of broad jets in identifying new physics signa- 
tures over QCD background. It is also disappointing be- 
cause MadGraph generated HELAS amplitudes [4] can be 
computed very fast on GPU (Graphic Processing Unit) [5J 

EI. 

Some of other simulation packages such as HELAC 
and Alpgen[Sl[7] employ QCD off-shell recursive relations 
to produce multi-parton scattering amplitudes efficiently. 
Successful computation of recursively evaluated QCD am- 
plitudes on GPU has been reported[5]. It is therefore in- 
teresting to examine the possibility of implementing re- 
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1 In the case of purely gluonic processes, even gg — > 5g pro- 
cess cannot be evaluated on MadGraph [3] . 



cursive relations for gluon currents in MadGraph without 
sacrificing its benefits. In this paper, we examine the use 
of the recursive amplitudes in the color-flow basis, since 
corresponding HELAS amplitude package has been devel- 
oped in ref . [5] ■ The main purpose of this paper is to show 
that this approach can be accommodated in MadGraph 
with explicit examples. 

The outline of this paper is as follows. In section [21 
we briefly review the color-flow decomposition and off- 
shell recursive relations. We then discuss the implemen- 
tation of those techniques in MadGraph, using its "user 
mode" option[3J in section[3] In section[3]and[5j we explain 
how we evaluate the color-summed amplitude squared and 
show the numerical results of n-jet production cross sec- 
tion. Section |5] is devoted to conclusions. 



2 The color-flow basis and off-shell recursive 
relations 

In this section we review the color-flow decomposition^ 
of QCD scattering amplitudes and the off-shell recursive 
relation [TU] of gluonic currents in the color flow basis. 



2.1 The color-flow decomposition 

First, we briefly review the color-flow decomposition. The 
Lagrangian of QCD can be expressed as 

C = -\{rn\i^u)) + '¥{i-fWi - m6}} <fc, (1) 
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where 



+ iJ^{(A^(A u )i - (A)j(AOfe} , (2) 



(3) 
(4) 
(5) 



Upper and lower indices are those of 3 and 3 representa- 
tions, respectively. Note that the gluon fields, (-4 P )* [i, j = 

1, 2, 3), are renormalized by a factor of y/2, and hence the 
coupling g is divided by v2- At this stage, not all the nine 
gluon fields are independent because of the traceslessness 
of the SU(3) generators, 



(AJl + (A,)l + (A,)l = 0. 
Here we introduce the "Abelian gluon" , 



(B 



IDj 



V2B L 



(6) 



(7) 



This is essentially the gauge boson , B^, of U(l) subgroup 
of U(3), combined with its generator, 5j/V2N, which is 
normalized to 1/2. We then rewrite eq. ([!} by adding and 
subtracting the Abelian gluon, such that the QCD La- 
grangian is expressed as 



V2 



V2 



(B, 



11)0 



V2 



- m5] 



with 



Mi = Wj + w 3 



(8) 



Here all the nine gluons, (G^)j, are now independent while 
the covariant tensors keep the same form 
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Fig. 1. QCD Feynman rules in the color-flow basis. The stan- 
dard Feynman diagrams are shown in the left and the color-flow 
diagrams are shown in the right. The wavy lines in the left and 
the dashed lines in the right denote the propagating Abelian 
gluon. The gluonic three point vertex has an additional non- 
cyclic permutation in color indices, while the four point vertex 
has additional five non-cyclic permutations. 



In this basis, these nine gluons have one-to-one correspon- 
dence to a set of indices of 3 and 3 representations of U (3), 
(G/j.)j, and we address them and this basis as 
U(3) gluons and the color-flow basis, respectively, in the 
following discussions. Although the definition of the co- 
variant tensor (flU)) contains self-coupling terms for U(3) 
gluons, the Abelian gluon contribution actually drops out 
in the sum. Accordingly, the Feynman rules for U(3) glu- 
ons give directly the SU (3) (QCD) amplitudes for purely 
gluonic processes [5]. 

We list the Feynman rules derived from the Lagrangian 
dSJ in Fig.[TJ It should be noted that both the propaga- 
tor and the quark vertex of the Abelian gluon have an 
extra —1/N factor because of the color-singlet projector, 



(9) Sj/VN. The negative sign of this factor directly comes 
from the sign of the Abelian gluon kinetic term in the 
Lagrangian. We can evaluate QCD amplitudes by using 
these Feynman rules. 

For a n-gluon scattering process, it is written as 

M(ng)= Y 5*/ M ■ ■ ■ 5 l ° {1) 

<J£S„/Z„ 

x A a (a{l),a{2),--- ,tr(n)), (12) 
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Fig. 2. All 6 color flows of the process GjlGj. 
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obtained from the corresponding n U(3) gluon scatter- 
ing process. A a (a(l),a(2), ■ ■ ■ , cx(n)J are gauge-invariant 

partial amplitudes, which depend only on particle mo- 
menta and helicities represented simply by the gluon in- 
dices, 1, • • • , n, and the summation is taken over all (n— 1)! 
non-cyclic permutations of gluons. A set of n Kronecker 
delta's in each term gives the "color flow" of the corre- 
sponding partial amplitudes. As an example, we list in 
Fig. [5] all six color flows indicated by Kronecker delta's for 
the gg -> gg process, -> Qf s Q)^ . 

On the other hand, if a scattering process involves 
quarks, contributions from the Abelian gluon do not de- 
couple. Therefore, we have to add all its contributions to 
the total amplitude. For processes with one quark line, 
such contributions come from Abelian gluon emitting dia- 
grams. The total amplitude for qq — > ng process becomes 

M(qq + ng) 



<res n 



ctGS„ v 7 

x ^(«,«7-(l)r-- ,<r(n~ o-(n)J 
' o-es„ v 7 

X Al(q,a(l),--- ,a(n-2),q; a(n-l),a(n)J 



■n) % 



■5l^Al{q,q; 



(13) 



where A k a [q, <r(l), • • • , <r(n — k), q; cr(n — fc + 1), • • • , cr(n) 

is the partial amplitude with k Abelian gluons, and q and 
q denote the momenta and the helicities of the quark and 
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Fig. 3. Color flows of qq — > ng process. For each color flow, 
one of the corresponding Feynman diagrams is shown in the 
parenthesis. The term "perms." stands for the other color flows 
obtained by permutations among gluons. 



the anti-quark, respectively. The summation is taken over 
all the n! permutations of n gluons. The first term in 
the r.h.s. of eq. (TT3"|) gives contributions from n E/(3) glu- 
ons, and the other terms give contributions from [n — k) 
U(3) gluons and k Abelian gluons, summed over k = 1 
to n, where (— 1/N) k comes from the —1/N factor in the 
Abelian gluon vertex as depicted in Fig. [TJ Corresponding 
color flows of those partial amplitudes are shown in Fig. [3] 

For processes with two or more quark lines, the Abelian 
gluon can be exchanged between them, giving the ampli- 
tudes with at least — 1/N suppression. Therefore, the to- 
tal amplitude of processes with two quark lines such as 
<7i(?2 — > <7i<72 + n gluons becomes 

M(qiq~[ + q 2 q2 + ng) 

n 

= EE 

er€S„ r=0 

\(tfii /cM") xM^+iA ( sfi2 *\ 

x Al(qi,<r(l),--- ,a(r),gS\qa,a(r + !),••• ,a(n),qi) 



— ( fi lq2 A 1 "'™ 3 x l '<x+i) 



jV(-r) jV(r-l) i<Ji 



X-B°Ui,o-(l),--- ,cr(r),gi | g 2 ,o-(r + 1),- • ■ ,o-(n),g 2 j 

+ (partial amplitudes with external Abelian gluons) . 

(14) 

A°(qi,(r(l),-~ ,a(r),qt\q 2 ,a(r+l),--- ,a(n),qi) (or B° a ) 

denotes a partial amplitude with n U(3) gluons, where two 
sets of arguments divided by "|" belong to two different 
color-flow chains; one starts from gi (qi in the initial state) 
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Fig. 4. Sample diagrams for the qiq2 — > qiqigg process. Each 
Feynman diagram in left-hand-side is accompanied by the cor- 
responding color-flow diagram in the right-hand-side. 



and ends with cfe (<?2 in the final state), and the other 
starts from (72 and ends with <7i , which are given explicitly 
by two sets of delta's before the partial amplitude. The 
difference between A° a and is that the former consists 
of U (3)-gluon-exchange diagrams while the latter consists 
of Abelian-gluon-exchange ones. Here we write down ex- 
plicitly the partial amplitudes with external U (3) gluons, 
A and B°, whereas those with external Abelian gluons 
should be added as in the previous one-quark-line case. For 
illustration, some typical diagrams for the q\q2 — > 91 5253 
process are shown in Fig. [4] 



2.2 Recursive relations of off-shell gluon currents 

Next, we mention the off-shell recursive relations in the 
color-flow basis. We define a (n + 1) -point off-shell gluon 
current J^(l, 2, • • • ,n; x) recursively as 

J**(Jfc; x) = 8f*8f k e' 1 (k) (1 < fc < n), (15) 



J M (/,/ + l,--- ,m; x) (l<l<m<n) 
—i 



m—l 



— ^ee v^(Pi, k ,p k+ x, m ,-Pi,, 



x J V {1, ■■■ , fc; y) J p {k + 1, • • • , m; z) 

111 — 2 m—l 

E E E stfz^vr»j^---My) 

y,z,w k=l q=k+l 



xJ p (k + l, ■■■ ,q; z)J a -(q + !,-■■ ,m;w)j, (16) 

where the numbers in the argument of currents represent 
gluons, and their order respects the color flow as in the 
partial amplitudes, A a , in eq. (|12[) ; the argument after 
semicolon denotes the off-shell gluon. As an initial value, 
the 2-point gluon current, J M (fc; x), is defined as the po- 
larization vector of the gluon k, e M (fc). Pij(i < j) is the 
partial sum of gluon momenta: 



Pi,j =Pi+ P1+1 + 



Vh 



(17) 



where pi to pj are defined as flowing-out momenta. V£ p>k 
and V 4 iyp<TM denote the three- and four-point gluon vertices 
from the Feynman rule of Fig.[T] 

V^{p, q, r )=-i^={{q- pyf + (r - g) V" 

+ (P-W} 1 (18) 

v u P an _ iL.(25* , V" - afg*" ~ g v>l 9 pc )- (19) 

The summation over y, z and w in eq. (|16[) are taken over 
the 3 and 3 indices, (j y ,i y ),(jz,iz) an d (jw,iw)i of off- 
shell gluons in currents J v , J p and J a , respectively. For 
illustration, the 5-point off-shell current is shown explicitly 
in Fig.EJ 

Helicity amplitudes of pure gluon processes in the color- 
flow basis are obtained from off-shell gluon currents as 

{n— 2 n-1 
EE E ttt v r p (Pi,k,p^Pi + ^ 
x,y,z k=l l=k+l 

x J M (1, • • • , fc; x)J u (k + 1, • • • ,1; y)J p {l + 1, • • • , n; z) 

n— 3 71— 2 n— 1 

E E E tttt^ 



E 

x,y,z,w k=l l=k+l m=l+l 



x J M (1, • • • , fc; x)J u {k + 1, • • • , Z; y)J p (l + 1, • • • , m; z) 
xJ CT (m + l,--- ,n;w) I. (20) 




Fig. 5. The 5-point gluon off-shell current. Numbers denote 
gluon momenta and helicities, and the curly line with x symbol 
denotes the off-shell gluon. 
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Fig. 6. An example of the 4-point off-shell current (the exter- 
nal leg with the cross symbol is off-shell), where an on-shell 
gluon in the current (labeled as 3) is replaced by an off-shell 
gluon attached to a quark line. 

Starting from the external wave functions (|15p , the recur- 
sive relation computes partial amplitudes very effectively 
by using the HELAS code. When the off-shell line "x" 
in Fig. [5] is set on-shell, the same set of diagrams give the 
complete helicity amplitudes in the color-flow basis. 

3 Implementation in MadGraph 

In this section, we discuss how we implement off-shell re- 
cursive relations in MadGraph in the color-flow basis and 
how we generate helicity amplitudes of QCD processes. 

3.1 Subroutines for off-shell recursive formulae 

First, we introduce new HELASfjE] subroutines in Mad- 
Graph which make n-point gluon off-shell currents and 
n-gluon amplitudes in the color-flow basis, according to 
eq. dTtj)) and eq. ([20)1 . respectively Although the expres- 
sion eq. (|16[) gives the off-shell gluon current made from 
n on-shell gluons, in HELAS amplitudes any input on- 
shell gluon wave functions can be replaced by arbitrary 
off-shell gluon currents with the same color quantum num- 
bers. Shown in Fig. [6] is an example of such diagrams that 
are calculated by the HELAS subroutine for the 4-point 
off-shell gluon current. Thanks to this property of the 
HELAS amplitudes, we need to introduce only two types 
of subroutines: one which computes helicity amplitudes of 
n-gluon processes via eq. (1201) and the other which com- 
putes off-shell gluon currents from n-point gluon vertices 
via eq. (|16p . We name the first type of subroutines as 
gluon# and the second type as jgluo#, where # denotes 
the number of external on-shell and off-shell gluons. 

The number of new subroutines we should add to Mad- 
Graph depends on the number of external partons (quarks 
and gluons) in QCD processes. Processes with n-partons 
can be classified as those with n gluons, those with (n— 2) 
gluons and one quark line, those with (n — 4) gluons and 
two quark lines, and so on. 

In the color-flow basis, the first class of processes with 
n external gluons are calculated by just one amplitude 
subroutine, gluon n. 

For the second class of processes with (n — 2) gluons 
and one quark line, we need up to (n — l)-point off-shell 
current subroutines, jgluo# with # = 4 to n — 1. This 



gluon# 


#-gluon amplitude in the color-flow basis 


jgluo# 


off-shefl gluon current from (#—1) external glu- 




ons in the color-flow basis 


ggggcf 


4-gluon amplitude from the contact 4-gluon ver- 




tex in the color-flow basis 


jgggcf 


off-shell gluon current from the contact 4-gluon 




vertex in the color-flow basis 


j ioaxx 


off-shell Abelian gluon current from a quark pair 




in the color-flow basis 



Table 1. New HELAS subroutines added into MadGraph. 
gluon* and jgluo# compute the #-gluon amplitude and the 
off-shell gluon current from the #-gluon vertex, respectively, 
by using the off-shell recursive formulae in the color-flow basis. 
We use gluon* with # = 4 to 7 and jgluo# with # = 4 to 6 
in this study. Two subroutines for the contact 4-gluon vertex 
(ggggcf and jgggcf ) are introduced to sum over the two chan- 
nels (s and t or s and u) for a given color flow, according to 
the Feynman rule in Fig.[T] We also add the off-shell Abelian 
gluon current subroutine from a quark pair (j ioaxx). 



is because the largest off-shell gluon current appears in 
the computation of diagrams where (n — 2) on-shell glu- 
ons are connected to the quark line through one off-shell 
gluon, which can be computed by the (n — l)-point off- 
shell current and the qqg amplitude subroutine, iovxxx. 
Note that the same diagram can be computed by the off- 
shell gluon current subroutine made by the quark pair, 
jioxxx, and the (n— l)-point gluon amplitude subroutine, 
gluon (n-1) . This type of redundancy in the computation 
of each Feynman diagram is inherent in the HELAS am- 
plitudes, which is often useful in testing the codes. 

For n-parton processes with (n — 4) gluons and two 
quark lines, we need up to (n — 2)-point off-shell current 
subroutines. By also introducing multiple gluon amplitude 
subroutines up to (n — 2)-point vertex, the maximum re- 
dundancy of HELAS amplitudes is achieved. Likewise, for 
n-parton processes with m quark lines, we need up to 
(n — m)-point off-shell current or amplitude subroutines. 

We list in Table [T] the new HELAS subroutines we 
introduce in this study. The subroutine gluon# evaluates 
a #-gluon amplitude in the color- flow basis, and jgluo# 
computes an off-shell current from (#— 1) external gluons. 
Since we consider up to seven parton processes, gg — > 5g, 
uu — > 5g and uu — ¥ uu3g, in this study, we use gluon4 to 
gluon7 and jgluo4 to jgluo6. 

In addition to these subroutines that computes am- 
plitudes and currents recursively, we also introduce three 
subroutines: ggggcf, jgggcf and j ioaxx. Two of them, 
ggggcf and jgggcf, evaluate an amplitude and an off- 
shell current from the contact 4-gluon vertex, following 
the Feynman rule of Fig.[TJ Although the amplitude sub- 
routine and the off-shell current subroutine for the 4-gluon 
vertex already exist in MadGraph, ggggxx and jgggxx, we 
should introduce the new ones which evaluate the sum of s- 
and t-type or s- and u-type vertices for a given color flow, 
since the default subroutines compute only one type at a 
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time, jioaxx computes an off-shell Abelian gluon current 
made by a quark pair and is essentially the same as the 
off-shell gluon current subroutine, jioxxx, in the HELAS 
library [J] except for an extra —N factor: 



jioaxx = —N x jioxxx. 



(21) 



Note that introducing this —N(= —1/N x N 2 ) factor is 
equivalent to summing up contributions from all Abelian 
gluon propagators, as we discussed in sectional We show 
the codes of ggggcf , jgggcf , gluon5 and jgluo5 in Appendix. 



3.2 Introduction of a new Model: CFQCD 

Next, we define a new Model with which MadGraph gen- 
erates HELAS amplitudes in the color-flow basis since the 
present MadGraph computes them in the color-ordered 
basis 11 . A Model is a set of definitions of particles, their 
interactions and couplings; there are about twenty pre- 
set Models in the present MadGraph package Ej such as 
the Standard Model and the Minimal SUSY Standard 
Model[T2]. There is also a template for further extension 
by users, User Mode (usrmod). Using this template, we 
can add new particles and their interactions to MadGraph. 
We make use of this utility and introduce a model which 
we call the CFQCD Model. 

In the CFQCD Model, we introduce gluons and quarks 
as new particles labeled by indices that dictate the color 
flow, such that the diagrams for partial amplitudes are 
generated according to the Feynman rules of Fig.Q] We 
need one index for quarks and two indices for gluons, such 
as itfc, Ttk and gji. The Abelian gluon, 30, does not have 
an index. The index labels all possible color flows, and it 
runs from 1 to m, where 

m = the number of gluons + the number of quark lines. 

This number m is the number of Kronecker's delta sym- 
bols which dictate the color flow. As an example, let us 
consider the purely gluonic process 

g(l)g(2)^g(3)g(4)g(5)g(6)g(7), (23) 

for which we need seven delta's to specify the color flow: 

^& (*»)5£ wSS • (24) 

Here the numbers in parentheses label gluons whereas the 
numbers in the sub-indices of Kronecker's delta's count 
the color-flow lines, 1 to m — 7 as depicted in Fig. [7] In 
CFQCD, we label gluons according to their flowing-out, 3, 
and flowing- in, 3, color-flow-line numbers, such that the 
partial amplitude with the color flow (|24p is generated as 
the amplitude for the process 

532(2)321(1) ->ffi 7 (7) 576(6)<?6 5 (5) .g54(4)<?43(3). (25) 



This is the description of the process (123)) in our CFQCD 
Model, and we let MadGraph generate the corresponding 
partial amplitudes, such as A(l, ■ ■ ■ , n) in eq. (|20p. as the 



9(1) 
i(D«(i) 



im 9P) 




Fig. 7. The color-flow diagram of the process Q23p with the 
color flow (|24[) . The number, 1 to 7, counts the Kronecker's 
delta's, or the color-flow lines. 

i(5)*(5) j(4)*(4) j(3)H3) J(2)«(2) j(l)«(l) 
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,,4 



- 3 



"2 W 



Fig. 8. The color-flow diagram for uu — > 5g process. The 
numbers in the diagram shows the color-flow-line number. 



helicity amplitudes for the process (|25l) . This index num- 
ber assignment has one-to-one correspondence with the 
color flow (pM|) . For instance, 532(2) in the process (|25|) 

denotes the contribution of the U(3) gluon G^ipi, ^2)^2) 

to the partial amplitude where its 3 index i(2) terminates 
the color-flow line 2, and the 3 index j(2) starts the new 
color-flow line 3. It should also be noted that we number 
the color-flow lines in the ascending order along the color 
flow, starting from the color-flow line 1 and ending with 
the color-flow line m. This numbering scheme plays an im- 
portant role in defining the interactions among CFQCD 
gluons and quarks. 

Let us now examine the case for one quark line. For 
the 5-jet production process 



ml ^ 5 (1) 5 (2) 3(3) 5 (4) 3(5), 



(26) 



the color-flow index should run from 1 to 6 according to 
the rule (l22l). Indeed the color flow 



corresponds to the process 



»(i) 



(27) 



Ui U 6 -t 565(5) 554(4) 343(3) 332(2) 321(1) 



(28) 

We show in Fig. [5] the color-flow diagram for this process. 
This is just that of the 6-gluon process cut at the 316 
gluon, where 316 gluon is replaced by the quark pair, u\ 
and 1*6. Shown in Fig. [H] are a few representative Feynman 
diagrams contributing to the process (|28D . Both the exter- 
nal and internal quarks and gluons in the CFQCD model 
are shown explicitly along the external and propagator 
lines. Following the MadGraph convention, we use flowing- 
out quantum numbers for gluons while the quantum num- 
bers are along the fermion number flow for quarks. Gluon 



Kaoru Hagiwara, Yoshitaro Takaesu: Generating QCD amplitudes in the color-flow basis with MadGraph 



7 



o I 

o 

o 
o 



P 
o 

o 

o 

o 





Fig. 9. A few representative Feynman diagrams in the color- 
flow basis for the color flow of Fig. [5] Arrows indicate the quan- 
tum number flows of gluons and quarks. 

J(5)«(5) j(4)?(4) j(3)«(3) j(2)i(2) j(l)Hl) 



ty ^ 



■Ju 



Fig. 10. A color-flow diagram for uu — > 5g process with 
Abelian gluons. The numbers in the diagram shows the color- 
flow-line number. 



propagators attached to quark lines are named by their 
color-flow quantum numbers along arrows. The diagram 
(a) contains only qqg vertices, (b) contains an off-shell 4- 
point gluon current, jgluo4, or a 4-point gluon amplitude, 
gluon4, and (c) contains either an off-shell 6-point gluon 
current, jgluo6, or a 6-point gluon amplitude, gluon6. 

In CFQCD, not only C7(3) gluons but also the Abelian 
gluon contributes to the processes with quark lines. For 
the uu —> 5g process ([26]) . 1 to 5 gluons can be Abelian 
gluons, go- If the number of Abelian gluons is k, the color 
flow reads 



(^•••OWr) 



*(5-fc + l) 

i(5-*+i) 



6-fe, 



(Si)}!'. (29) 



When k = 5, all the five gluons are Abelian, and the 
first (the right-most) color flow should be (^i)j", just as 
in uu — > 5 photons. In Fig. lIOl we show the color-flow 
diagram for the process (I2B")) with the color flow (|2^|) for 
k — 3 as an example. 

For the process with two quark lines, 



The CFQCD model computes the partial amplitude for 
the above color flow as the helicity amplitude for the pro- 
cess 

u x d 5 7/5^4343(3)332(2)^21(1) ■ (32) 

The other possible color-flow processes without Abelian 
gluons are 

uidi u 5 d 3 . 932(2)321(1). 954(3), (33) 
u x d 3 -> u 5 d 2 .921(1) .954(3) .943(2), (34) 
uid 2 -> u 5 di .954(3) 343(2) .932(1), (35) 

for k — 3, 2, 1, respectively. As in the single quark line case, 
all the external gluons can also be Abelian gluons, and 
we should sum over contributions from external Abelian 
gluons. For instance, if the gluon 3(3) in the process 
is Abelian, the color flow becomes 



(36) 



and the corresponding partial amplitude is calculated for 
the process 



uid 4 U4,d 3 332(2) 321(1)30(3) 



(37) 



in CFQCD. 

When there are more than one quark line, the Abelian 
gluon can be exchanged between two quark lines, and the 
color flow along each quark line is disconnected. For in- 
stance, the color flow 



H2) (x \Hi) 



(38) 



is obtained when the Abelian gluon is exchanged between 
the w-quark and d-quark lines. The corresponding partial 
amplitude is obtained for the CFQCD process 



uid 4 -> u 3 d 5 332(2) 321(1) 354(3). 



(39) 



Note that in the process (|39|) each flow starts from and 
ends with a quark pair which belongs to the same fermion 
line. 

We have shown so far that we can generate diagrams 
with definite color flow by assigning integer labels to quarks, 
9^, and gluons, gki, such that the labels k and I count the 
color-flow lines whose maximum number m is the sum of 
the number of external gluons and the numbers of quark 
lines (quark pairs), see eq. ([212]) . 



ud^udg(l)g(2)g(3), (30) 

the color- flow index should run from 1 to 5. All the pos- 
sible color flows are obtained from the comb-like diagram 
for the one-quark-line process shown in Fig. [31 by cutting 
one of the gluons into a quark pair, such as gk+i,k to dk+i 
and dk- Then the first color flow starting from ui ends at 
dk, and the new color flow starts with dk+i which ends at 
u m . An example of such color flow for k — 4 and m = 5 
reads 



3.3 New particles and their interactions in the CFQCD 
Model 

In Table [2[ we list all the new particles in the CFQCD 
model for m = 6 in eq. (f22|) . The list is shown in the 
format of particles . dat in the usrmod of MadGraph[2]. 
As explained above, the C/(3) gluons have two color-flow 
indices, gu, while the Abelian gluon, 30, has no color- flow 
index. They are vector bosons (type=V) , and we use curly 
lines (line = C) for U (3) gluons while wavy lines (line — W) 
for the Abelian gluon. The 'color' column of the list is used 
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particle 


ant i- p ar t ic le 


type 


line 


mass 


width 


color 


laoci 


code 


U(3) gluons 


















gl3 


g!3 


V 


C 


ZERO 


ZERO 


S 


g!3 


21 


gl4 


gl4 


V 


C 


ZERO 


ZERO 


s 


gl4 


21 


gl5 


gl5 


V 


C 


ZERO 


ZERO 


s 


gl5 


21 


gl6 


gl6 


V 


C 


ZERO 


ZERO 


s 


gl6 


21 


gai 


g21 


V 


C 


ZERO 


ZERO 


s 


g21 


21 


p-24 


s-24 
o 


v 


c 


ZERO 


ZERO 


s 


p-24 


21 


g 25 


g25 


V 


c 


ZERO 


ZERO 


s 


g 25 


21 


g26 


g26 


V 


c 


ZERO 


ZERO 


s 


g26 


21 


gb2 


go 2 


T J 

V 


ft 




Z.r.KU 


b 


gb2 


z 1 


g63 


g63 


V 


c 


ZERO 


ZERO 


s 


g63 


21 


g64 


g64 


V 


c 


ZERO 


ZERO 


s 


g64 


21 


g65 


g65 


V 


c 


ZERO 


ZERO 


s 


g65 


21 


Abelian gluon 


















go 


go 


V 


w 


ZERO 


ZERO 


s 


go 


21 


Quarks 


















ul 


ul~ 


F 


s 


ZERO 


ZERO 


s 


ul 


2 


u2 


u2~ 


F 


s 


ZERO 


ZERO 


s 


u2 


2 


u6 


u6~ 


F 


s 


ZERO 


ZERO 


s 


u6 


2 



Table 2. New particles of the CFQCD model when m = the number of gluons + the number of quark lines = 6. The list is 
shown in the format ofparticles.dat in the usrmod of MadGraph[2]. 



by MadGraph to perform color summation by using the 
color-ordered basis. Since we sum over the color degrees 
of freedom by summing the 3 and 3 indices (j's and i's) 
over all possible color flows explicitly we declare all our 
new particles as singlets (color = S]3- The last column 
gives the PDG code for the particles, and all our gluons, 
including the Abelian gluon, are given the number 21. 

All gluons, not only the Abelian gluon but also C/(3) 
gluons, are declared as Majorana particles (particle and 
anti-particle are the same) in CFQCD. We adopt this 
assignment in order to avoid generating gluon propaga- 
tors between multi-gluon vertices. Such a propagator is 
made from a particle coming from one of the vertices and 
its anti-particle from the other. Since we define the anti- 
particle of a U (3) gluon as the gluon itself, the color flow 
of the propagator should be flipped as shown in Fig. ??, 
according to the gluon naming scheme explained in the 
previous subsection. Therefore, CFQCD does not give di- 
agrams with gluon propagator in "color-flow conserving" 
amplitudes. For example, to form a c/13 gluon propaga- 
tor between two multi-gluon vertices, we need a 313 gluon 
coming out from one of the vertices and a (731 gluon from 
the other, according to the gluon naming scheme explained 
in the previous subsection. However, because a 1731 gluon 
is not an anti-particle of g\3 but a different particle in our 




2 If we declare all CFQCD gluons as octet and quarks 
as triplets, MadGraph generates color-factor matrices in the 
color-ordered basis for each process, which is not only useless 
in CFQCD but also consumes memory space. 



CF conserving 
qqg vertex 

gluon propagator 

CF flipping 
qqg vertex 

Fig. 11. 'Majorana' (7(3) gluon propagator in CFQCD. k and 
I denote the color-flow-line number. Since MadGraph propa- 
gator connects a flowing-out particle with a flowing-out anti- 
particle, it necessarily flips the color flow, as shown along the 
vertical lines. It either breaks the color-flow conservation, (a), 
and hence cannot appear inside pure gluon amplitudes, or it 
contributes to amplitudes with quark lines with the help of 
the color-flow flipping qqg vertex, denoted by a cross symbol 
in. (b). 



particle difinition, Table [2], 1713 gluon propagator can not 
be formed. U(3) gluon propagators attached to quark lines 
are allowed by introducing the color-flow non-conserving 
qqg couplings that effectively recover a color-flow connec- 
tion. See discussions on qqg vertices at the end of this 
section. 

In Tabled we list u-quarks, Uk, and its anti-particles, 
Uk, with the color-flow index k = 1 to 6. They are all 
fcmions (type = F) for which we use solid lines (line = S) 
in the diagram. Their colors are declared as singlets (color 
= S) as explained above. The list should be extended to 
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the other five quarks, such as dk and dk for down and 
ant i- down quarks. 

Before closing the explanation of new particles in CFQCD 
let us note that the number of U(3) gluons needed for m- 
gluon processes is m(m — 2). This follows from our color- 
flow-line numbering scheme, which counts the successive 
color-flow lines in the ascending order along the color flow 
as depicted in Fig.[7J This is necessary to avoid double 
counting of the same color flow. According to this rule, the 
only gluonic vertex possible for to = 3 is the one among 
3i3> 332 and 1721 ■ Although 1731 can appear for processes 
with to > 4, (712 and (723 can never appear. Generalizing 
this rule, we find that gki's with I = k + 1 mod(m) as well 
as I = k cannot appear and hence we need only to(to — 2) 
gluons in CFQCD. 

In this paper, we report results up to 5-jet production 
processes. Purely gluonic gg — > bg process has seven ex- 
ternal gluons, and to = 7 is necessary only for this process. 
According to the above rule, there is only one interaction 
for this process, which is the one among 1717, g^, 365, 354) 
343; 332 and 321- Therefore, we should add 1717 and gre to 
Table O in order to compute gg — » bg amplitudes. 

Next, we list all the interactions among CFQCD par- 
ticles. Once the interactions among particles of a user 
defined model are given, MadGraph generates Feynman 
diagrams for an arbitrary process and the correspond- 
ing HELAS amplitude code. In CFQCD, we introduce 
n-point gluon interactions and let MadGraph generate 
a code which calls just one HELAS subroutine (either 
gluon# or jgluo# with # = n, which makes use of the 
recursion relations of eqs. (f2"0"l) or (fTB)) . respectively) for 
each n-point vertex. In addition, we have quark-quark- 
gluon vertices. We list the interactions in the descending 
order of the number of participating particles in Table [3] 
tod 

First, we show the 7-point interaction in Table [3] in the 
format of interactions .dat of usrmod in MadGraph[2 . 
As discussed above, it is needed only for generating gg — ► 
bg amplitudes. The vertex is proportional to g^, the fifth 
power of the strong coupling constant, and we give the 
five couplings, cpll to cpl5, as 

G2 = -%, (40) 

according to the Feynman rules of Fig. [iflj whose types, 
typel to type5, are all QCD. 

The 6-point interactions appear in gg — ¥ Ag and also in 
qq — > bg process in this study. Again, only one interaction 
is possible among the six gluons, gki, with fc = 1 to 6 and 
I = k — 1 , as shown in Table 01 The coupling order is g s 
and the four couplings, cpll to cpl4, are G2 all with the 
type QCD. 

The 5-point gluon vertices appear in gg — > 3g, qq — ► 
43, qq — > bg and qq —> qq3g processes. The first two and 
the fourth processes have m = 5, for which the 5-point 
gluon vertex is unique as shown in the first row of Ta- 
ble El The process qq — > bg gives m = 6, and gluons 

3 The couplings in the HELAS codes are defined to be —i 
times the couplings of the standard Feynman rule. 



Table 3. 7-point vertices 



particles cpll ■ ■ ■ cpl5 typcl ■ 


■ ■ type5 


gl7 g76 g65 g54 g43 g32 g21 G2 ■ • ■ G2 QCD • 


■ • QCD 


Table 4. 6-point vertices 


particles cpll ■ ■ ■ cpl4 typcl ■ ■ ■ 


typc4 


gl6 g65 g54 g43 g32 g21 G2 ■ • ■ G2 QCD ■ ■ ■ 


QCD 



Table 5. 5-point vertices 



particles cpll cpl2 cpl3 typcl typc2 typc3 



gl5 g54 g43 g32 g21 G2 G2 G2 QCD QCD QCD 

g26 g65 g54 g43 g32 G2 G2 G2 QCD QCD QCD 



Table 6. 4-point vertices 



particles 


cpll 


cpl2 


typel typc2 


gl4 g43 g32 g21 


G2 


G2 


QCD 


QCD 


g25 g54 g43 g32 


G2 


G2 


QCD 


QCD 


g36 g65 g54 g43 


G2 


G2 


QCD 


QCD 


gl5 g53 g32 g21 


G2 


G2 


QCD 


QCD 


gl5 g54 g42 g21 


G2 


G2 


QCD 


QCD 


gl5 g54 g43 g31 


G2 


G2 


QCD 


QCD 



Table 3 to 6: G2= —g s /y/2 is the coupling of gluon vertices 
and defined as a real, according to the HELAS convention. 
The list is shown in the format of interactions.dat in the 
usrmod of MadGraph [2] . 



with the sixth color-flow line can contribute to the 5-point 
gluon vertex. Because of the ascending color-flow number- 
ing scheme, only one additional combination appears as 
shown in the second row of Table [5] These vertices have 
<7g order and we have cpl/s = G2 and typefc = QCD for 
k= 1,2,3. 

The 4-point gluon vertices appear in gg — ¥ gg (to = 4) , 
qq — > (m— l)g with m — 4 to 6, and in qq — > qq+ (m — 2)g 
with m — 4 and 5 in this study. As above, we have only one 
4-point gluon vertex for processes with m = A, which is 
shown in the first row of Table [S] For the process qq — > Ag 
(to = 5), one additional vertex appears as shown in the 
second row. In case of qq — > bg (to = 6), the ordering 
3 — > A — » 5 — > 6 also appears, and it is given in the third 
row. 

So far, we obtain multiple gluon vertices from the color- 
flow lines of consecutive numbers, corresponding to the 
color flows such as 

y + l-*y + 2->-"--*y + n, 

(0<y<m-n) (41) 
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Table 7. 3-point vertices 



particles 


cpl 


type 


q-q-q vertices 

j <j <3 






gl3 g32 g21 


G2 


QCD 


g24 g43 g32 


G2 


QCD 


g35 g54 g43 


G2 


QCD 


g46 g65 g54 


G2 


QCD 


gl4 k42 k21 

o o o 


G2 


QCD 


gl4 k43 e31 

O o o 


G2 


QCD 


gl5 g52 g21 


G2 


QCD 


gl5 g54 g41 


G2 


QCD 


p-oc; p-c;^ p-^o 


G2 


QCD 


g25 g54 g42 


G2 


QCD 


q-q-g vertices 






ul u3 gl3 


GG2 


QCD 


ul u3 g31 


GG2 


QCD 


u3 ul g31 


GG2 


QCD 


ul u4 gl4 


GG2 


QCD 


ul u4 g41 


GG2 


QCD 


u4 ul g41 


GG2 


QCD 


u6 u7 g76 


GG2 


QCD 


u7 u6 g76 


GG2 


QCD 


q-q-go vertices 






ul ul gO 


GGO 


QCD 


u2 u2 gO 


GGO 


QCD 


u6 u6 gO 


GGO 


QCD 



G2= -g s /V2, GG2= {-g s /V2, -g s /V2) and 
GG0= (g s /(\^2N), g a / (\/2N)) are the couplings of each vertex. 
G2 is defined as a real and GG2 and GGO are defined as two 
dimensional complex arrays, according to the HELAS 
convention. The list is shown in the format of 
interactions.dat in the usrmod of MadGraph [2]. 



for n-point gluon vertices. When there are two or more 
quark lines in the process, we can also have color flows 
which skips color-flow-line numbers, such as 

y + 1 — y + n — j/ + n + <f, 
(0 < y < m — n , — d), 
(l<d<m-n) (42) 

for n-point gluon vertices, where d counts the number of 
skips. Because gluon propagators do not attach to gluon 
vertices in CFQCD, this skip can appear only when two 
or more gluons from the same vertex are connected to 
quark lines. For example, in the fourth row of Table [51 the 
gluon <7i5 couples to the quark line, u§ — > u±, and then 
the gluon 353 couples to the other quark line, d^ — > ds, 
in the process ud — ¥ udgg. Likewise, 342 and (731 in the 



fifth and the sixth row, respectively, couple to the second 
quark line, d — > d, of the process. 

As for the 3-gluon vertices, ten vertices listed in Table[7] 
appear in our study. The first four vertices with successive 
color-flow numbers appear in processes with one quark 
line, qq — > (m — l)g with m = 3 to 6, and those with 
two quark lines, qq — > qq (m — 2)g with m = 4 to 5. The 
fifth and the sixth vertices starting with gl4 appears for 
qq — » qq (m — 2)g with m = 4, for which only one unit of 
skip (d = 1) appears, and also with m = 5. The last four 
vertices appear only in the m — 5 process with two quark 
lines, for which d — 2 is possible. In fact, the two vertices 
starting with gl5 contain g52 or g41 with d — 2 skips in 
the color-flow number. This completes all the gluon self- 
interactions in CFQCD up to 5-jet production processes. 

There are two types of qqg vertices in CFQCD: the cou- 
plings of U (3) gluons, gki, and those of the Abelian gluon, 
go- All the qqg couplings for the M-quark, u% to ue, are 
listed in Table [7] In the HELAS convention, the fermion- 
fermion-boson couplings are two dimensional complex ar- 
rays, where the first and the second components are the 
couplings of the left- and the right-hand chirality of the 
flowing-in fermion. According to the Feynman rules of 
FigfTJ the couplings are 



(43) 



for U(3) gluons and 



GGO = 



(44) 



for the Abelian gluon; Note the 
gluon coupling. 

U(3) gluons have interactions as 

ui iik gu, 



( g s g s 
KVZN' s/2N , 

1/N factor for the Abelian 



(45) 



where the fermion number flows from ui to Uk by emitting 
the out-going g^i gluon. All the diagrams with on-shell 
U(3) gluons attached to a quark line are obtained from 
the vertices (|45j) . Likewise, the Abelian gluon couples to 
quarks as 

u k u k g . (46) 

In CFQCD, we generate an U (3) gluon propagator be- 
tween a quark line and a gluon vertex and between two 
quark lines by introducing color-flow flipping vertices 



Uk ui g k i 



(47) 



as we discussed above. Here, the gluon gki is emitted from 
the quark line Uk — > ui, such that U(3) gluons can prop- 
agate between a quark line and a gluonic vertex and also 
between two quark lines when the other side of the qqg 
vertex is the color- flow conserving one, (|4"5)) . In order to 
exchange gluons between an arbitrary gluonic vertex and 
a quark line, these color-flow-flipped qqg vertices should 
exist for all U(3) gluons. However, since all U(3) gluons 
have the color-flow conserving vertices (|4"51 as well, we find 
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double counting of amplitudes where an U (3) gluon is ex- 
changed between two quark lines. This double counting 
can be avoided simply by discarding color-flow conserv- 
ing qqg vertices (|45[) for k < I as shown in Table 7. 

Now we exhaust all the interactions needed in the cal- 
culation of partial amplitudes and ready to generate di- 
agrams and evaluate total amplitudes in the color-flow 
basis. 



4 Total amplitudes and the color summation 

In this section, we discuss how we evaluate the total am- 
plitudes, eqs. (H2), CCU and (JUJ), with the CFQCD model 
and how we perform the color summation of the total am- 
plitude squared. 

4.1 gg — > ng processes 

We consider pure gluon processes first. As discussed in 
section [21 the total amplitude of a n-gluon process is ex- 
pressed as eq. (|12p . and it is the sum of several partial 
amplitudes. Because color factors, Kronecker's delta's, for 
partial amplitudes is either 1 or 0, the total amplitude for 
a given color configuration (a color assignment for each 
external gluons) consists of a subset of all the partial am- 
plitudes in eq. ([12")) . Therefore, in order to evaluate the 
total amplitude for a given color configuration, we should 
find all possible color flows for the configuration, compute 
their partial amplitudes and simply sum them up. 
As an example, let us consider a n = 5 case, 

g(l)g(2) -> 5(3)5(4)5(5). (48) 

The color configuration is expressed by sets of 3 and 3 
indices, j and i, respectively, for each gluon; 

(jl,h)(l),{h,k)(2),--- ,(j5,«5)(5)- (49) 

Here the subscripts of each index pair label the gluon g(k), 
k = 1 to 5, with the momentum pk and the helicity 
For instance, let us examine a color configuration 

(1,1) (1) ,(1,1) (2) ,(2,1) (3) ,(3,2) ( 4),(1,3) (5 ). (50) 

One of the possible color flows that gives the above color 
configuration is 

(*«)isw5S!w5S(*)5S(*o;S' ( 51 ) 

whose associated amplitude can be calculated as the am- 
plitude for the CFQCD process 

ff32(2)s2i(l)-> 515(5)554(4)543(3). (52) 

MadGraph generates the corresponding HELAS ampli- 
tude code, which calls the 5-gluon amplitude subroutine, 
gluon5. 

There is also another color flow 



for the color configuration ([5H)1 . The corresponding partial 
amplitude can be obtained by one of the (n — 1)! permu- 
tations of (n — 1) gluon momenta and helicitics: 

{^(5),e^(4),^(3),^(2),^(l)} 
->{^(2),^(5),^(4),^(3),^(1)} ) (54) 

where e M (i) denotes the wave function of the external 
gluon g(i). 

When all the partial amplitudes for each color flow 
are evaluated, we sum them up and obtain the color- 
fixed total amplitude, eq. (fl"2")). for this color configura- 
tion. Therefore, for pure gluon process, we generate the 
HELAS amplitude once for all to evaluate the color-fixed 
total amplitude. The total amplitude is then squared and 
the summation over all color configurations is done by the 
Monte-Carlo method. 



4.2 qq — > ng processes 

Next, we discuss the processes with one quark line, qq —> 
ng. The procedure to compute the color-summed ampli- 
tude squared is the same, but we should take into ac- 
count the Abelian gluon contributions for processes with 
quarks. As shown in Fig.[3J the Abelian gluons appear as 
an isolated color index pair which have the same number. 
Therefore, we can take into account its contribution by re- 
garding the independent gluon index pair as the Abelian 
gluon. 

As an example, let us consider the process 

uu^ 5(1)5(2)5(3)5(4) (55) 

with the color configuration: 

(l,l) u ,(l,l) (1) ,(2,l) (2) ,(3,2) (3) ,(l,3) (4) , (56) 

where the parenthesis with the subscript 'it' gives the color 
charges of the u-quark pair; (j u ,i u ) u denotes the annihila- 
tion of the u-quark with the 3 charge, j u , and the u-quark 
with the 3 charge, i u . Although this color configuration 
is essentially the same as ([ST)]) , and hence has the color 
flows like (|5Tj) and ([53)1 . we have an additional color flow 
for this process: 

(*«);s^&^)3S(*«)5S^)S J B (57) 

Here the index pair of the first gluon, (1, l)m in (|56")h 
forms an independent color flow, (^Tjny which corre- 
sponds to the Abelian gluon, 5o(l)- The CFQCD process 
for this color flow reads 

u{u 4 -> 543(4) 532(3) 521(2) 50(1). 

4 As a convention, we always count the color-flow-lines from 
the one which starts from the 3 index of the quark, j q , for 
processes with quark lines. 
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Summing up, there are three color flows for the color con- 
figuration (|56l) : 



(a) 



(6) 



(a) ; u — y 1 — y 2 — y 3 — y 4 — y u, 

(b) :u-»2->-3->4->-l->-7l, 

(c) :u^2^3^4^u, 1 



1. 



(58a) 
(58b) 
(58c) 



The partial amplitude A\, for the color flow (b) is obtained 
from that of the color flow (a), A a , by a permutation of 
gluon wave functions as in the case of gg 3g amplitudes 
for the color configuration (fBTJf . The partial amplitude for 
the color flow (c) is calculated with the Feynman rules of 
Fig-Ill an d the total amplitude M. is obtained as 



M = A a 



(59) 



where the —1/N factor comes from the Abelian gluon as 
shown in eq. ([TB")) . When more than one gluon have the 
same color indices, jk = ik in (jk,ik)(k), more Abelian 
gluons can contribute to the amplitude, and the factor 
(— 1/N) m appears for the partial amplitude with m Abelian 
gluons. In the CFQCD model introduced in section [31 
those factors are automatically taken into account in the 
HELAS code generated by MadGraph. 



4.3 qq -> qq(n — 2)g processes 

Finally, let us discuss the processes with two quark lines. 
As shown in Fig.|H there are two independent color flows 
because there are two sets of quark color indices, (j qi , i qi ) 
and (jq 2 ,iq 2 )- Nevertheless, we can show color configura- 
tions and find possible color flows in the same way as in 
the qq — ¥ ng cases. Let us consider the process 



ud -> ud g(l)g(2)g(3) 
with the color configuration 



(60) 



(l,l) u ,(l,l) d ,(2,l) (1) ,(3,2) (2) ,(l,3) (3) (61) 

for illustration. Color charges of u-quarks are in the paren- 
thesis {j u , i u )u, and those of d-quarks are in (jd, id)d as in 
the qq — > ng processes. All the possible color flows are 



(u^l^2^3^ d,d^u) 

cb):(*«)5- w («.);s( tf »)5S(*»)5!, l) (*oS: 
w («*)3" w (*»)5S ^fi' 

(u — > 1 — > 2 — > 3 — > u,d — > d) 
(u — > u,d — > f ^ 2 ^ 3 — > d) 



(62a) 
(62b) 
(62c) 
(62d) 



We show in Fig.rT2] one representative Feynman diagram 
for each color flow, (a) to (d). The amplitudes for each 



| 
p 

I 
« p 



IS 



(d) 



<# 

4> 4> 4> 

3 2 1 

Fig. 12. Feynman diagrams for the color flow (a) to (d). We 
show one of the Feynman diagrams for each color flow. 



color flow, A a to Ad, are calculated as amplitudes for the 
corresponding CFQCD processes: 

(a) : u\ d 5 -> u 5 d 4 543(3) 532(2) 521(f), (63a) 

(b) :uid 2 ^ u 5 di 554(3) 543(2) 532(f), (63b) 

(c) :uid 5 -> u 4 d 5 543(3) 532(2) 521(f), (63c) 

(d) :uid 2 ^ ui d 5 554(3) 543(2) 532(f)- (63d) 

For the CFQCD processes (c) and (d), the color flow lines 
starting with the u- and d-quarks terminate at the same 
quarks. Such amplitudes are generated by an exchange 
of the Abelian gluon, as shown by the representative di- 
agrams in Fig. Q21 The total amplitude is now obtained 
as 

M = A a +A b -—A c -—A d , (64) 



N 



N 



according to eq. (I14|) . The color summation of the squared 
amplitudes \Ai\ 2 is performed by the MC method just as 
in the pure gluon case. 



5 Sample results 

In this section, we present numerical results for several 
multi-jet production processes as a demonstration of our 
CFQCD model on MadGraph. We compute n-jet produc- 
tion cross sections from 55 — > ng, uu — > ng and uu — > 
uu -f- (n — 2)5 subprocesses up to n — 5 in the pp colli- 
sion at y/s — 14 TeV. We define the final state cuts, the 
QCD coupling constant and the parton distribution func- 
tion exactly as the same as those of ref. [3], so that we 
can compare our results against those presented in ref. [3] , 
which have been tested by a few independent calculations. 
Specifically, we select jets (partons) that satisfy 



1 7ft I < 2.5, 

pr(j), PT ]k > 20 GeV, 



(65a) 
(65b) 
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gg -> ng 














uu - 

( 


> uu +(n-2)g 
v/o Abelian) 






(1 Abelian) _ 




( 


uu -> ng - 
w/o Abelian) 









2 3 4 



The number of jets 

Fig. 13. Total cross sections of gg — > ng (upper), uu — > uu + 
(n — 2)g (middle) and uu — > ng (lower) in fb scale for pp 
collisions at yfs = 14 TeV as given in Table [8] Results shown 
by dashed lines and dotted lines for uu and uu subprocesses 
are obtained when Abelian gluon contributions are ignored and 
one Abelian gluon contributions are included, respectively. 



more than 30% for n = 2, about 13% for n = 3, while 
about 10% for n —A and 5. 

The bottom lines show results for the subprocess uu — > 
ng. As above, the solid line gives the exact cross sections, 
while the dashed line and the dotted line give the results 
when contributions from Abelian gluons are ignored and 
one Abelian gluon contributions are included, respectively. 
Unlike the case for uu ^ uu+ (n — 2)g subprocesses, the 
Abelian gluon contributions remain at 30% level even for 
n = 5. 

Before closing this section, we would like to give two 
technical remarks on our implementation of CFQCD on 
MadGraph. First, since the present MadGraph[2 does 
not allow vertices among more than 4 particles, we add 
HELAS codes for 5, 6, and 7 gluon vertices by hand to 
complete the MadGraph generated codes. This restric- 
tion will disappear once the new version of MadGraph, 
MGE0, is available, since MG5 accepts vertices with arbi- 
trary number of particles. Second, we do not expect dif- 
ficulty in running CFQCD codes on GPU, since all the 
codes we developed (see Appendix) follow the standard 
HELAS subroutine rules. 



where r\j and pt{j) are the pseudo-rapidity and the trans- 
verse momentum of the parton-j, and PT h is the smaller of 
the relative transverse momentum between parton-j and 
parton-A:. We use CTEQ6L1 parton distribution functions [13 
at the factorization scale Q = 20 GeV and the QCD cou- 
pling a s (Q = 20 GeV)j^g = 0.171. Phase space integra- 
tion and summation over color and helicity are performed 
by an adaptive Monte Carlo (MC) integration program 

BASES ng. 

Results are shown in Table [5] and FigfTHl 
In Table [51 the first row for each n-jet process gives 
the exact result for the n-jet production cross section, 
while the second row shows the cross section when we ig- 
nore all the Abelian gluon contributions. The third row 
shows the results where we include up to one Abelian 
gluon contributions: one Abelian gluon emission ampli- 
tudes without an Abelian gluon exchange for uu — > ng 
and uu — y uu -)- (n — 2)g processes, and Abelian gluon ex- 
change amplitudes without an Abelian gluon emission for 
uu — ^ uu + (n — 2)g processes. All the numerical results for 
the exact cross sections in Table [5] agree with those pre- 
sented in ref. [3] within the accuracy of MC integration. In 
Fig. 1131 the multi-jet production cross sections are shown 
for n = 2,3,4, and 5 jets in units of fb. The upper line 
gives the results for the subprocess gg — > ng. The middle 
lines show those for the subprocess uu — >• uu + (n — 2)g; 
their amplitudes are obtained from those of the subprocess 
ud —> ud+ (n — 2)g outlined in this study, simply by anti- 
symmetrizing the amplitudes with respect to the two ex- 
ternal quark wave functions. The solid line gives the exact 
results, while the dashed line gives the results when all the 
Abelian gluon contributions are ignored. The dotted line 
shows the results which include up to one Abelian gluon 
contributions, although it is hard to distinguish from the 
solid line in the figure. Despite their order 1/N C suppres- 
sion, contributions of Abelian gluons can be significant; 



6 Conclusions 

In this paper, we have implemented off-shell recursive for- 
mulae for gluon currents in the color-flow basis in Mad- 
Graph and have shown that it is possible to generate QCD 
amplitudes in the color-flow basis by introducing a new 
model, CFQCD, in which quarks and gluons are labeled 
by color-flow numbers. We have 

- introduced new subroutines for off-shell recursive formu- 
lae for gluon currents, the contact 3- and 4-point gluon 
vertices in the color flow basis and the off-shell Abelian 
gluon current, 

- defined new MadGraph model: the CFQCD Model, 

- generated HELAS amplitudes for given color flows and 
calculated the color-summed total amplitude squared, and 

- showed the numerical results for n-jet production cross 
sections (n < 5). 

Although we have studied only up to 5-jet production 
processes in this paper, it is straightforward to extend the 
method to higher n-jet production processes. 
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No. of jets 


gg -¥ ng 


uii — > ng 




uu — > uu (n 


-2) S 


2 


(3.19 ± 0.00) x 10 11 


(2.90 ± 0.00) 


X 


10 7 


(2.67 ± 0.01) 


x 10 8 






(4.23 ±0.01) 


x 


10 7 


(3.65 ±0.00) 


x 10 8 






(2.97 ±0.00) 


X 


10 7 


(2.67 ±0.01) 


x 10 8 


3 


(2.61 ± 0.00) x 10 10 


(1.84 ±0.00) 


X 


10 6 


(5.88 ± 0.03) 


x 10 7 






(2.65 ±0.01) 


X 


10 6 


(6.66 ± 0.04) 


x 10 7 






(1.96 ±0.01) 


X 


10 6 


(5.93 ±0.03) 


x 10 7 


4 


(5.81 ± 0.01) x 10 9 


(4.52 ±0.10) 


X 


10 5 


(2.78 ±0.02) 


x 10 7 






(6.17 ±0.01) 


X 


10 5 


(3.03 ±0.03) 


x 10 7 






(4.78 ±0.01) 


X 


10 5 


(2.80 ±0.02) 


x 10 7 


5 


(4.57 ±0.01) x 10 9 


(1.57 ±0.03) 


X 


10 5 


(1.54 ±0.01) 


x 10 7 






(2.13 ±0.01) 


X 


10 5 


(1.73 ±0.01) 


x 10 7 






(1.57 ±0.04) 


X 


10 s 


(1.56 ±0.01) 


x 10 7 



Table 8. Total cross sections of gg — > ng, uu — > ng and uu — > uu + (n — 2)g (n < 5) in fb scale for pp collisions at yfs = 14 TeV, 
when jets satisfy 1 77^,- 1 < 2.5, pr(j) > 20 GeV and pr jk > 20 GeV for the smaller of the relative transverse momentum between 
two jets. Results in the second row and the third row of each n-jet cross section are obtained when we ignore Abelian gluon 
contributions and include one Abelian gluon contributions, respectively. Abelian gluons do not contribute to purely gluonic 
processes. 



Appendix: Sample codes for off-shell currents and 
amplitudes 

In this Appendix, we list HELAS codes for the con- 
tact 4-point gluon vertex subroutines, ggggcf and jgggcf , 
which sums over contributions with definite color-flow. In 
addition, we list HELAS codes for the 5-gluon amplitude 
subroutine, gluon5, and the 5-point off-shell gluon current 
subroutine, jgluo5, as examples of the recursive multi- 
gluon vertices introduced in section [3l 



A1 - gggg cf 

subroutine ggggcf (gl ,g2 ,g3 ,g4 ,g, vertex) 
implicit none 

complex*16 gl(6) ,g2(6) ,g3(6) ,g4(6) .vertex 
complex*16 dvl(0:3) ,dv2(0:3) ,dv3(0:3) ,dv4(0:3) 
complex*16 dvertx , vl2 , vl3 , vl4 , v23 , v24 , v34 
real*8 g 

dvl(0)=dcmplx(gl(D) 
dvl(l)=dcmplx(gl(2)) 
dvl(2)=dcmplx(gl(3)) 
dvl(3)=dcmplx(gl(4)) 
dv2(0)=dcmplx(g2(D) 
dv2(l)=dcmplx(g2(2)) 
dv2(2)=dcmplx(g2(3)) 
dv2 (3) =dcmplx (g2 (4) ) 
dv3(0)=dcmplx(g3(D) 
dv3(l)=dcmplx(g3(2)) 



dv3 (2) =dcmplx (g3 (3) ) 
dv3 (3) =dcmplx (g3 (4) ) 
dv4 (0) =dcmplx (g4 ( 1 ) ) 
dv4(l)=dcmplx(g4(2)) 
dv4 (2) =dcmplx (g4 (3) ) 
dv4 (3) =dcmplx (g4 (4) ) 

vl2= dvl(0)*dv2(0)-dvl(l)*dv2(l)-dvl(2)*dv2(2) 
& -dvl(3)*dv2(3) 

vl3= dvl(0)*dv3(0)-dvl(l)*dv3(l)-dvl(2)*dv3(2) 
& -dvl(3)*dv3(3) 

vl4= dvl(0)*dv4(0)-dvl(l)*dv4(l)-dvl(2)*dv4(2) 
& -dvl(3)*dv4(3) 

v23= dv2 (0) *dv3 (0) -dv2 ( 1) *dv3 ( 1 ) -dv2 (2) *dv3 (2) 
& -dv2(3)*dv3(3) 

v24= dv2 ( ) * dv4 ( ) - dv2 ( 1 ) *dv4 ( 1 ) -dv2 ( 2 ) *dv4 ( 2 ) 
& -dv2(3)*dv4(3) 

v34= dv3 ( ) * dv4 ( ) - dv3 ( 1 ) *dv4 ( 1 ) -dv3 ( 2 ) *dv4 ( 2 ) 
& -dv3(3)*dv4(3) 

dvertx = (-vl4*v23+2 . d0*vl3*v24-vl2*v34) 

vertex = dcmplx (dvertx) *(g*g) 

return 
end 

A2 - jggg cf 
subroutine jgggcf (wl ,w2,w3,g, jggg) 

implicit none 



double complex wl(6) ,w2(6) ,w3(6) , jggg(6) 
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double complex dwl(0:3) ,dw2(0:3) ,dw3(0:3) 
double complex j j (0 : 3) ,dv,w32 ,wl3,w21 ,dg2 
double precision q(0:3,)q2,g 

double precision rOne 
parameter ( rOne = l.OdO ) 

jggg(5) = wl(5)+w2(5)+w3(5) 
jggg(6) = wl(6)+w2(6)+w3(6) 



dwl(0) 




dcmplx 


(wi(D) 


dwl(l) 




dcmplx 


(wl(2)) 


dwl(2) 




dcmplx 


(wl(3)) 


dwl(3) 




dcmplx 


(wl(4)) 


dw2(0) 




dcmplx 


(w2(l)) 


dw2(l) 




dcmplx 


(w2(2)) 


dw2(2) 




dcmplx 


(w2(3)) 


dw2(3) 




dcmplx 


(w2(4)) 


dw3(0) 




dcmplx 


(w3(D) 


dw3(l) 




dcmplx 


(w3(2)) 


dw3(2) 




dcmplx 


(w3(3)) 


dw3(3) 




dcmplx 


(w3(4)) 



q(0) = -dble(jggg(5)) 
q(l) = -dble(jggg(6)) 
q(2) = -dimag(jggg(6)) 
q(3) = -dimag(jggg(5)) 

q2 = q(0)**2-(q(l)**2+q(2)**2+q(3)**2) 
d g 2 = g*g 

dv = r0ne/dcmplx(q2) 

w32 = dw3(0)*dw2(0)-dw3(l)*dw2(l) 
& -dw3(2)*dw2(2)-dw3(3)*dw2(3) 

wl3 = dwl(0)*dw3(0)-dwl(l)*dw3(l) 
& -dwl(2)*dw3(2)-dwl(3)*dw3(3) 

w21 = dw2(0)*dwl(0)-dw2(l)*dwl(l) 
& -dw2(2)*dwl(2)-dw2(3)*dwl(3) 

jj(0) = dg2*(-dwl(0)*w32+2d0*dw2(0)*wl3 

& -dw3(0)*w21) 

jj(l) = dg2*(-dwl(l)*w32+2d0*dw2(l)*wl3 

& -dw3(l)*w21) 

jj(2) = dg2*(-dwl(2)*w32+2d0*dw2(2)*wl3 

& -dw3(2)*w21) 

jj(3) = dg2*(-dwl(3)*w32+2d0*dw2(3)*wl3 

& -dw3(3)*w21) 

jggg(l) = dcmplx (j j (0) *dv) 
jggg(2) = dcmplx(j j (1) *dv) 
jggg(3) = dcmplx(jj(2)*dv) 
jggg(4) = dcmplx(j j (3)*dv) 



15 

return 
end 

A3. gluon5 

subroutine gluon5 (wl,w2,w3,w4,w5,g, vertex) 
implicit none 
integer i 

complex*16 w(6,12,12) ,wl(6) ,w2(6) ,w3(6) ,w4(6) 
complex*16 w5(6) ,z(55) ,wx(6,55) , vertex 
real*8 g 

vertex=(0d0,0d0) 



do i=l,6 




w(i,l,l) 


= wl(i) 


w(i,2,2) 


= w2(i) 


w(i,3,3) 


= w3(i) 


w(i,4,4) 


= w4(i) 


w(i,5,5) 


= w5(i) 


enddo 





do i=l,3 

call jggxxx(w(l,i,i) ,w(l , i+1 , i+1) ,g, 
& w(l,i,i+D) 
enddo 

do i=l,2 

call jggxxx(w(l,i,i+l) ,w(l , i+2 , i+2) ,g, 
& wx(l,l)) 

call jggxxx(w(l,i,i) ,w(l , i+1 , i+2) ,g,wx(l,2)) 

call jgggcf (w(l,i,i) ,w(l,i+l,i+l) , 
& w(l,i+2,i+2) ,g,wx(l,3)) 

call sumw(wx,3,w(l,i,i+2)) 
enddo 

call gggxxx(w( 1,1,3) ,w(l,4,4) ,w(l,5,5) ,g,z(l)) 
call gggxxx(w( 1,1,2) ,w(l,3,4) ,w(l,5,5) ,g,z(2)) 
call gggxxx(w(l,l,l) ,w(l,2,4) ,w(l,5,5) ,g,z(3)) 
call ggggcf (w( 1,1, 2) ,w(l,3,3) ,w(l,4,4) , 

& w(l,5,5) ,g,z(4)) 

call ggggcf (w(l, 1,1) ,w(l,2,3) ,w(l,4,4) , 

& w(l,5,5) ,g,z(5)) 

call ggggcf (w(l,l,l),w(l,2,2),w(l,3,4), 

k w(l,5,5),g,z(6)) 
do i=l,6 

vertex = vertex+z(i) 
enddo 



return 
end 
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A4. jgluo5 

subroutine jgluo5(wl,w2,w3,w4,g, jgluon5) 
implicit none 



do i=l,3 

call jggxxx(w(l,i,i) ,w(l , i+1 , i+1) ,g, 
k w(l,i,i+l)) 
enddo 

do i=l,2 

call jggxxx(w(l,i,i+l) ,w(l , i+2 , i+2) ,g, 

& wx(l,l)) 

call jggxxx(w(l,i,i) ,w(l , i+1 , i+2) ,g,wx(l,2)) 
call jgggcf (w(l,i,i) ,w(l,i+l,i+l) , 

& w(l,i+2,i+2) ,g,wx(l,3)) 

call sumw(wx,3,w(l,i,i+2)) 
enddo 

call jggxxx(w(l,l,3) ,w(l,4,4) ,g,wx(l,l)) 
call jggxxx(w(l,l,2) ,w(l,3,4) ,g,wx(l,2)) 
call jggxxx(w(l,l,l) ,w(l,2,4) ,g,wx(l,3)) 
call jgggcf (w(l, 1,2) ,w(l,3,3) ,w(l,4,4),g, 

& wx(l,4)) 

call jgggcf (w(l, 1,1) ,w(l,2,3) ,w(l,4,4),g, 

& wx(l,5)) 

call jgggcf (w(l, 1,1) ,w(l,2,2) ,w(l,3,4),g, 

& wx(l,6)) 
call sumw(wx,6, jgluon5) 

return 
end 
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integer l 

complex*16 w(6,12,12) ,wl(6) ,w2(6) ,w3(6) 
complex*16 w4(6) ,wx(6,55) , jgluon5(6) 
real*8 g 

do i=l,6 

jgluon5(i) = (0d0,0d0) 

w(i,l,l) = wl(i) 

w(i,2,2) = w2(i) 

w(i,3,3) = w3(i) 

w(i,4,4) = w4(i) 
enddo 
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